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The numerical solution for the B = 2 static soliton of the SU{2) Skyrme model shows a profile 
function dependence which is not exactly radial. We propose to quantify this with the introduction 
of an axially symmetric oblate ansatz parametrized by a scale factor d. We then obtain a relatively 
deformed bound soliton configuration with Mb^2/Mb^i = 1.958. This is the first step towards to 
description of B > 1 quantized states such as the deuteron with a non-rigid oblate ansatz where 
deformations due to centrifugal effects are expected to be more important. 
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II. INTRODUCTION 

The Skyrme model I*] is a nonlinear effective field theory of weakly coupled pions in which baryons emerge as 
localized finite energy soliton solutions. The stability of such solitons is guaranteed by the existence of a conserved 
topological charge interpreted as the quantum baryon number B. More specifically, Skyrmions consist in static pion 
field configurations which minimize the energy functional of the Skyrme model in a given nontrivial topological sector. 
The model is partly motivated by the large- A^c QCD analysis there are reasons to believe that once properly 

quantized, a refined version of the model could accurately depict nucleons as well as heavier atomic nuclei with mesonic 
degrees of freedom 0,0,0 in the low energy limit. 

In the lowest nontrivial topological sector B — 1, the Skyrmion is described by the spherically symmetric hedgehog 
ansatz which reproduces experimental data with an accuracy of 30% or better 0. However, this relative success 
radically contrasts with the situation encountered in the B > 2 sector, where the hedgehog ansatz is not the lowest 
energy configuration and would not give rise to bound state configurations Moreover, pioneering numerical 

investigations of Verbaarschot ^3 clearly indicate that the B = 2 Skyrmion is not spherically symmetric, but rather 
possesses an axial symmetry reflected in its doughnut-like baryon density. Further inspection of this numerical 
solution, in particular the profile function, suggests that the classical biskyrmion may be represented by an oblate 
field configuration. Yet, most of the trial functions used to describe such a solution assume a decoupling from the 
angular degrees of freedom, i.e. F{r,6,(f>) = F{r), as this is the case for the instanton-inspired ansatz proposed in 
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or in the early variational approach |l2l Il3l | for example. On the other hand, some axially symmetric solutions 
were analysed in the B — 1 sector, by |14| and respectively, to include possible deformations due to centrifugal 
effects undergone by the rotating Skyrmion and account for the quadrupole deformations of baryons. Here, our aim 
is to extend the work on oblate Skyrmions in TFl to describe dibaryon states. The classical static oblate solution 
introduced in this manner will provide a quantitative estimate of the axial deformation, which is different from a 
uniform scaling in a given direction as performed in |14| . It should also provide an adequate ground to perform the 
quantization of the B = 2 soliton. 

There has been several attempts to decribed the angular dependence of the B > I solutions which is much more 
complicated than the hedgehog form in i? — 1. Fortunately, a few years ago, Houghton et al. came up with an 
interesting ansatz based on rational maps. The rational map ansatz provides a simple alternative compared to the full 
numerical study of the angular dependence of the baryon density distribution of multiskyrmions. It also yields static 
energy predictions in good agreement with numerical solutions for several values of B. The most interesting feature 
of this method remains without doubt that fundamental symmetries of multiskyrmions can easily be implemented 
in the ansatz solutions. This provides a clever way to identify the symmetries of the exact solutions, which are not 
always apparent, and in some cases, adequate initial solutions for lengthy numerical calculations. Close as it may be, 
the rational map ansatz remains an approximation and in some cases, more accurate angular ansatz have been found. 
For example, Houghton and Krusch |l6| slightly improved the mass approximation of the biskyrmion by relaxing 
the requirement of holomorphicity imposed on rational maps. However, the profile function defined in this work still 
solely depends on r. As may seems evident, some accuracy still may be gained by introducing a more appropriate 
parametrization of the soliton shape function. Recently, loannidou et al. |l7l | obtained similar results by introducing 
an improved harmonic maps ansatz where the profile function depends on radial and polar degrees of freedom as well. 
However, they had to deal with a complicated second-order partial differential equation. 

From these considerations, we propose a B = 2 oblate solution based on rational maps, which could be understood 
as the rational maps solution proposed in 0,^^ with the radial dependence F{r) replaced by an oblate form /(r/). 
Consequently, the soliton undergoes a smoothly flattening along a given axis of symmetry. The parameter d provides a 
measure of the scale at wich the deformation becomes important while the solution preserves the angular dependence 
given by the rational maps scheme for the B — 2 case. This choice is obviously consistent with the toroidal baryon 
density of the B = 2 Skyrmion. Implementing the oblate ansatz to the model, we first inte grat e analytically the 
angular degrees of freedom. This explains why other angular ansatz such as those in and |17| were not chosen; 
they led to complications. The second step involve solving the remaining nonlinear ordinary second-order differential 
equation resulting from the minimization of the static energy functional with respect to the profile function f{r]). 
Thereby, the parameter d is set as to minimize the static energy, i.e. the mass of the soliton. Although the method 
applies to higher baryon numbers and other Skyrme model extensions, the analysis is restricted here to i? = 2 
Skyrmion for SU (2) Skyrme model. 

In the next section, we present the axially symmetric oblate ansatz for the SU(2) Skyrme model introducing the 
oblate spheroidal coordinates. In section III, we briefly describe rational maps and show how they can be used in 
the context of static oblate biskyrmions. A discussion of the numerical results follows in the last section, where 
we also draw concluding remarks about how the oblate-like solution could be a good starting point to perform the 
quantization of the non-rigid B = 2 soliton as the deuteron. 



III. STATIC OBLATE MULTISKYRMIONS 



Let us first introduce the oblate spheroidal coordinates (77, 0, (j)), which are related to Cartesian coordinates through 



{x,y,z) = d(cosh77sin0cos(/), coshTysin^sini/), sinhT^cos^), (1) 

so a surface of constant rj correspond to a sphere of radius d cosh 77 flattened in the 2-direction by a factor of tanhyy. 
For small 77, these surfaces are quite similar to that of pancakes of radius d whereas when 77 is large, they become 
spherical shells of radius given by {d/2)e^K Note that taking the double limit d 0, 77 — > cx3 such that r always 
remains finite, one recovers the usual spherical coordinates. Thus, the choice of the parameter d establishes the scale 
at which the oblateness becomes significant. Finally, the element of volume reads 

dV = -d^ cosh 77 (sinh^ 77 + cos^ 9) drj d{cos 9) dcj). (2) 
Neglecting the pion mass term, the chirally invariant Lagrangian of the SU{2) Skyrme model just reads 

C = -^TtL^L^ + ^Tt[L^,L,]' (3) 
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where = Wd^U with U € SU{2). Here, is the pion decay constant and e is sometimes referred to as the 
Skyrme parameter. In order to implement an oblate solution, let us now replace the hedgehog ansatz 

U ^e'^ F(r) (4) 

by the static oblate solution defined as follow 

U = e" ■ ^ ^("^ (5) 

where the Tk stand for the Pauli matrices while rj is the standard unit vector rj = V77/IV77I. More explicitly, this unit 
vector is simply 

^ (coshr/ sin 8 cos $, cosh?7 sin sin(f>, sinhr/ cos 0) 

V= , ^ = , (6) 

J cosh rj — cos^ 0) 

As will become apparent in the next section, we consider the case where Q = 0(9) and <!> = <f>(0), i.e. 8 and $ 
depend only on the polar angle 6 and the azimuthal angle (f> respectively. Furthermore, /(ry), which determines the 
global shape of the soliton, plays the role of the so-called profile function. In that respect, the oblate ansatz is clearly 
different from a scale transformation along one of the axis [lj| • As in its original hedgehog form, the field configuration 
U constitutes a map from the physical space TZ'^ onto the Lie group manifold of SU (2). Finite energy solutions require 
that this SU{2) valued field goes to the trivial vacuum for asymptotically large distances, that is U{rj 00) l2x2- 
The expression for the static energy density is 

£^£2+£i = -^Tr{L,L,) + ^Tr L,f) (7) 
so, after substituting the oblate ansatz, the mass functional can be written as 



M[/(77), 8(0), $(</>)] ^ j dV {M2 + Mi) (8) 

with 

Al2-^(|Vr7p/'2 + sinVi^) (9) 

and 



Mi = ^ f 2| V?7|2 /'2 k + sin^ / (k'' - kahKab) ) ■ (10) 
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Here, the notation is lighten by the use of the Kab matrix defined as: 

Kab = ^iVa^iVb (11) 



K = Tr{Kab) = Kaa. (12) 

Introducing an auxiliary variable for convenience, 

S = cos^ 8 sin^ 8 + 8'^ cosh^ 77 sinh^ r] (13) 

one easily deduces that 

k=(^^ + (sin- 8/ sin- ^) \ ^^^^ 

\ (cosh 77 — cos^ 8)2 d2(cosh ry — cos^ 8) / 

k^ K k 2 .,2 |Vr/|-I](sin-8/sin-g) 
a (cosh 77 - cos2 8)3 

However, before minimizing the mass functional with respect to the chiral angle /(jy), in view to get the static 
configuration of the soliton, one must specify an angular dependence in Q{9) and $(0). This is the subjet of the next 
section, after a brief recall of some basic features related to the rational maps ansatz. 



and 
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IV. OBLATE SKYRMIONS AND RATIONAL MAPS 

Formally, a rational map of order N consists in a 5^ — > 5^ liolomorphic map of the form 

R.i.) = ^ (16) 

where p and q appear as polynomials of degree at most N. Moreover, these maps are built in such a way that p or 
q is precisely of degree N . It is also assumed that p and q do not share any common factor. Any point z on 
is identified via stereographic projection, defined through z = tan(|) e*"^. Thus, the image of a rational map R{z) 
applied on a point z of a Riemann sphere corresponds to the unit vector 

^R^Y+Jlh^ (25Ri?Jv(^),23i?jvW,l-|i?ivP) (17) 

which also belongs to a Riemann sphere. The link between static soliton chiral fields and rational maps ^ follows 
from the ansatz 

[/j^(r,z) = e"«-^^('') (18) 

inasmuch as F{r) acts as radial chiral angle function. To be well defined at the origin and at r — > oo, the boundary 
conditions must be i^(0) = kir where k is an integer and F{oo) = 0. The baryon number is given hy B = Nk where 
N — max(degp, degg) is the degree of Rn{z). We consider only the case fc = 1 here, so B — N. 

By analogy with the nonlinear theory of elasticity !T3| , Manton has showed that the static energy of Skyrmions 
could be understood as the local stretching induced by the map U : Ti? . In this real rubber-sheet geometry, the 

Jacobian J of the transformation provides a basic measure of the local distortion caused by the map U . This enables 
us to build a symmetric positive definite strain tensor defined at every point of Ti? as 

D.j^XjJ ^-^TriUL,). (19) 

This strain tensor £>, changing into DO under orthogonal transformations, comes with three invariants expressed 
in terms of its eigenvalues A^, A| and A|: 

Tr D = \l + \l + \l (20) 

\{Tr Df - \Tr = A?A^ + AiA^ + A^A^ (21) 

detD = \l\l\l. (22) 

Since it is assumed that geometrical distorsion is unaffected by rotations of the coordinates frame in both space and 
isospace, the energy density should remain invariant and could be written as a function of the basic invariants as 
follows 

8 = a{\\ + \l + \l)+(3 + A^A^ + A^A^) (23) 
where a and j3 are parameter depending on and e, while the baryon density is associated with the quantity 

- (24) 
27r^ 

In this picture, radial strains are orthogonal to angular ones. Moreover, owing to the conformal aspect of i?jv(z), 
angular strains are isotropic. Thereby, it is customary to identify 

Ai = ~F'{r) (25) 

and 

sini^(r) l + |zp 



A2 — A3 — 



1 + \Rn{z) 



di?Ar(z) 



dz 



(26) 
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Thus, substituting these eigenvalues in (|23() and integrating over physical space yields 



Mn^ I drl^£i'+-£i' ] (27) 



with 



=(f'' + 2N^^^] (28) 



7N _ lnAri^'2^''P - Sin^^ 
^2 



£^ = ( 2NF"—^ +1^—— ) (29) 



and 
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If 2idzdz f 1 + \z\^ 

In = 



An J |l + |z|2|2 Vl + |i?Ar(2)| 



dz 



(30) 



wherein ^"^"^^^J-i^i corresponds to the usual area element on a 2-sphere, that is sin^ d9 d(f>. At first glance, one sees 
that radial and angular contributions to the static energy are clearly singled out. 
Now, focussing on the N = 2 case, the most general rational map reads 



Mz) = — J— — r- (31) 

llZ^ + UZ + A 



However, imposing the exact 2— torus symmetries (axial symmetry and rotations of 180° around Cartesian axes), as 
expected from numerical analysis [13 , restricts the general form above to this one 

R.iz) = (32) 

It has been showed in (5(| that T2, ^^nd thus the mass functional alike, exhibits a minimum for a = 0. Then, one must 
conclude that the most adequate choice of a rational map for the description of the biskyrmion solution boils down 
to R2{z) = z"^ . Recasting this map in term of angular variables (0, 0) (8, $), we get 

z R^iz) = tan(® ) e'* = tan^i^) e^"^ 



e ^Q{9)^ arcsin 
(j) $(0) = 2(j) 



sin^ I 



1 + cos2 



It is easy to verify that for B = 1, the rational map is simply R{z) — z, or Q{9) — 9 and $(0) = 4>, and one recovers 
the energy density in 

We shall assume here that this angular function still holds in the oblate picture in H9ll5fl . After analytical angular 
integrations are performed, the mass of the oblate biskyrmion can be cast in the form 

47re f"" /J 1 \ 

M=— dTi[ E2 + -~eA (33) 



A Jo \2 " 4d ^ 
with 

^2 =X2i/'^ + Al22sinV (34) 

and 

£4 = X4i/'2sinV + X42sinV- (35) 

The explicit expressions of the density functions Mij{iii), which are reported in Appendix A, follow from straightfor- 
ward but tedious calculations. Adopting the same conventions as in jlSj , for the sake of comparison, we also rescaled 
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the deformation parameter as d = eFTrd/2^/2, with e = l/-\/2e and A = The values of F-^ = 129 MeV and 

e = 5.45 are chosen to coincide with those of 0- Now, the chiral angle function f{r]) can be determined from the 
minimization of the above functional, i.e. requiring that SM[f{ri)] = . Thus static field configuration must obey the 
following nonlinear second-order ODE: 

= /" (2dcoshr]+—.M4ism^n + ( ^7^41 sin / cos / 
+ /' ^2dsinh77 + ^-^41 sin^ / 



— (17^22 sin / cos / — ;Al42 sin"^ / cos / 
d 



(36) 



Here, the primes merely denote derivatives with respect to rj. Solving numerically for several values of d, we obtain 
the set of chiral angle functions of Figure 1 . When d is small, we recover exactly the solution of the N = 2 rational 
map ansatz. Let us stress that increasing d enforces a continuous displacement of the function /(rj) which induces a 
smooth deformation of the soliton. 




FIG. 1: Oblate biskyrmion chiral angle for several values of d. From left to right, we have d = 0.5, 0.2,0.1, 0.05, 0.01, 0.005 
and 0.001. 
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FIG. 2: The mass of the oblate biskyrmion solution as a function of d. The mass reaches its minimal value for d = 0.225. 



In Figure 2, we plot the mass of the oblate biskyrmion as a function of d. The mass of the biskyrmion passes 
trough a minimum for a finite non-zero value of the parameter d. This is a clear indication that the oblate solution is 
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energically favored. Note again that in the limit d ^ 0, we reproduce the mass value found in |5|,ll3 with the rational 
maps ansatz and profile function with radial dependence F(r). 



Numerical calculations carried out some years ago by Verbaarschot [lOj and almost concurrently by Kopeliovich and 
Stern establish that in the Skyrme model the mass ratio R2/1 = Mb=2/Mb=i is 1.92. The B = 2 oblate solution 
also represents a bound state of two solitons since its mass is lower than twice the mass found in the B = 1 sector 
From our calculations, we get that the mass of the static oblate biskyrmion, being minimized for d ~ 0.225, is 
16.066^ or 1689.5 MeV. The parameter d — 0.225 provides a measure of how the the B — 2 solution is flattened. 
Hence, the mass ratio of our flattened solution turns out to be R2/1 = 1.958. Comparing to other ansatz for _B = 2 
solutions, it is fairly smaller than that predicted by the familiar hedgehog ansatz with boundary conditions F{0) = 2t: 
and F(oo) — 0, since then R2/1 > 3 |3] or and still better than the hedgehog-like solution with (j) ^ 2(j) proposed 
as in ;9.j|, whose mass ratio is -R2/1 = 2.14 . Note that none of these solutions are stable solutions since R2/1 > 2. 
Let us mention that Kurihara et al. achieved a mass ratio of 1.94 using a different angular parametrization. 
However, there are no obvious physical grounds for their angular trial function and it remains that rational maps 
are far more superior when it comes to depict the symmetries of the B > 1 solutions. Our results still represent a 
slight improvement over those obtained in the original framework of rational maps, i.e. R2/1 — 1.962 where the 
chiral angle is strictly radial. Hence the oblate solution which depends of a spheroidal oblate coordinate rj, captures 
more exactly the profile shape of the biskyrmion than the original rational maps ansatz does although both rely on 
rational maps. The relatively small improvement also suggests that a better ansatz for classical B = 2 static solution 
would require a different choice for the angular dependence. In that regard, 13 and [13] both achieved a mass ratio 
of 1.933 by dropping the constraint on the rational maps to be holomorphic. These alternatives remain very difficult 
to implement for an oblate field configuration. 

It is worth emphasizing that the procedure presented here generalizes to any baryon number and any choice of 
angular ansatz consistent with multiskyrmions symmetries although analytical angular integration may become more 
cumbersome if not impossible in those cases. Similarly, the approach can be generalized to other Skyrme-like effective 
Lagrangians. 

In this paper we only investigated the classical B = 2 static solution but in principle, the full solution requires a 
quantization treatment to account for the quantum properties of dibaryons. The most standard procedure consists in 
a semiclassical quantization using collective variables. It only adds simple kinetic terms to the Hamiltonian but these 
energy contributions should partially fill the energy gap between the {I — Q, J — 1) deuteron mass (1876 MeV) and 
that of our deformed B — 2 static Skyrmion (1689.5 MeV). So, even if our analytical oblate ansatz is not necessarily 
the lowest static energy solution for the B = 2 Skyrmion, the optimization of the oblateness parameter d should 
prove adequate to take into account the soliton deformation due to centrifugal effects, as for the B = 1 case [Tsj . 
Thus, following such a procedure, we can expect that the properly quantized biskyrmion solution would provide a 
good starting point for the description of the low energy phenomenology of the deuteron 0, 121I |2^ . The problem of 
quantization of the oblate biskyrmion solution is an important topic in itself and will be addressed elsewhere. 
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Performing angular integrations in the oblate static energy functional, we get the following Mij {rj) density functions: 
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Mil = 



L{r]) (2 cosh"'^^ T] — 15 cosh"'^° 77 + 45 cosh^ 77 — 86 cosh^ rj + 124 cosh'* 77 — 92 cosh^ 77 - 



TT 



- 53 cosh^ 



2(cosh^ ?7 



- 4 cosh'' 77 + 8 cosh^ 77 - 4)2 
120 cosh^ 77 — 136 cosh'' 7/ + 76 cosh^ 77 — 16) 



(8 cosh*" 77 - 53 cosh° 7; ~r u 

2\/ cosh^ 7/ — 1 (cosh^ 7/ — 4 cosh* 7; + 8 cosh^ ~ 4)' 

(cosh^ 77 



8 cosh^ rj arctan ( , * 



5 cosh* 77 + 7 cosh^ '7 — 3) 



•\/ cosh^ 77 — 1 (cosh® 7; — 4 cosh'' 77 + 8 cosh^ 7; — 4)^ 

(2 cosh** 7; — 11 cosh^ 7; + 30 cosh^ 7; — 40 cosh'"^ rj + 28 cosh^ 7; — 8 cosh 7;) 
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(2 cosh* 7/ + cosh^ ?y — 7) £(77) 



2 cosh7;(2 cosh* 77 — 11 cosh^ ^7 + 7) 
(cosh^ r\—l) 



2 cosh^ 7; ycosl? 
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